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Abstract 

A method for fast simulation of the Cherenkov Hght generated by electromagnetic 
showers is described. The parametrization for the longitudinal profile is used and 
fluctuations and correlations of the parameters are taken into account in a consistent 
way. Our method dramatically reduces the CPU time and its results are in rather 
good agreement with a full Monte Carlo simulation. 
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1 Introduction 



The flux of tlic Clicrenkov pliotons from a sliower is typically simulated by 
tracking all secondary particles of the shower down to the Cherenkov thresh- 
old energy Ec- For the electrons in water £'c=0.78 MeV. The computer time 
needed for simulations of this type increases linearily with shower energy. 
Therefore a large amount of CPU time is required for simulation of a high en- 
ergy shower. However large samples of events are necessary to obtain results 
with a sufficient statistics. The parametrization of the longitudinal Cherenkov 
light profile of the showers is one of the methods to speed up simulation. 

A simple algorithm for parameterized showers has been successfully used 
for the simulation of the UAl calorimeter [1] . The simulation of the longitudi- 
nal energy profile of electromagnetic showers was based on fitting the param- 
eters of a Gamma distribution to the average shower profile. Later the shape 
fluctuations of individual showers were systematically taken into account for 
the simulation of calorimeter built for the HI experiment at HERA [2]. 

For a fast generation of the Cherenkov light from showers with energies 
Eq >10 GeV we have developed the program FLG. The program SIMEX [3] 
has been used for the full Monte Carlo (MC) simulation of the electromag- 
netic showers and computation of the Cherenkov light emitted by the shower 
particles. 



2 Parametrization of electromagnetic showier 

It is well known that the average longitudinal profile of elecromagnetic 
showers is reasonably well described by the Gamma distribution 
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where t is a shower depth in units of radiation length, 
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The parameters /3 and {a — l)/t are coefficients of absorbtion and creation of 
the shower particles, respectively. There is a relation between the location of 
the shower maximum t^ax. and the parameters a and /3, 

W = (" - (4) 



which follows from Eq.(l). The energy dependence of tmax can be approxi- 
mated by the expression 

imax^lnf— )+C, (5) 



where C—-0.5 for e-induced shower, C—0.5 for 7- induced shower and £=75.5 
MeV is the critical energy for water. 

The longitudinal shower profile f{t) reveals a useful scaling property when 
the depth of a shower is expressed in units t^iax as follows: T—t/tmax- Then 
Eq.(l) can be rewritten as 

dr i [a) 



and P' = t,„s.^f3. 

The total distance traveled by all charged shower particles L (track length 
in the units of radiation length) is proportinal to the shower energy Eq. 

L oc Eo/e. (7) 



3 Parametrization of the Cherenkov light profile of shower 



As we deal with underwater/underice detector measuring the Cherenkov 
radiation we will consider the average longitudinal Cherenkov light distri- 
bution emitted by the shower particle only. Total number of the Cherenkov 
photons from the shower, N^^^, is proportional to the track lenght L. But ac- 
cording to Eq.(7) the track lenght is proportional to the shower energy, and 
we can write 

NZ, = D,Eo. (8) 



So, the average longitudinal Cherenkov light distribution fchi^) can be parametrized 
as follows 

/ch(T) - D,f{T), (9) 
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where 

oo 

J U{T)dT = Nl, = D^Eo. (10) 



Then the number of the Cherenkov photons emitted by the shower particles 
in the layer At at the depth r is 

N(r, At) = ^ [7(0., (3'(t + Ar)) - ^{a, (3'r)] , (11) 



where 

X 

-f{a,x) = j r-^exp{-t)dt (12) 




is the incomplete Gamma function. 

A full MC simulation of the electromagnetic shower in water ( the radiation 
lenght Xo=36.1 cm) has been done Tising the program SIMEX. The value for 
the tracking energy threshold was set to 1 MeV. We obtain that the number 
of photons in the wavelenght 300 nm < A < 600 nm per MeV of the shower 
energy is 

= NljEo ^ 163 photons/MeV. (13) 



4 Simulation of the Cherenkov light profile taking into account the 
fluctuations of showers 



Realistic simulation, however, requires the simulation of individual show- 
ers. Fluctuations of the parameters in Eq.(9) obtained from the averaged 
Cherenkov light profile does not necessarly lead to a correct description of 
the fluctuations of individual showers. Assuming that also individual shower 
profiles can be approximated by the Gamma distribution 

. E.,(Airp^, (M) 

r(«j) 



the fluctuations and correlations can be taken into account consistently. The 
index i indicates that the function describes an individual shower i with the 
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parameters ctj and The ctj and P'^ can be calculated from the first fj and sec- 
ond af moments of the distribution Eq.(14) for each single SIMEX-simulated 
shower: 

C(i = {n/aif and (5[ = fi/a1. (15) 



The statistics of the full MC simulation was 1000 events for each fixed 
energy Eq in the region 10 GeV < Eq < 10"^ GeV and 200 events for ^0=10^ 
GeV. The parameters ctj and (3'^ are normal-distributed and such that the 
means (oii), (PI) and their fiuctuations ctq and ap can be determined and 
parametrized as a function of the shower energy. The average values (aj) and 
{f3^) vs shower energy are shown in Fig.l. They can be approximated by a 
logarithmic energy dependence. The correlation of the ctj and is given by 



and slowly increases with shower energy from p ~ 0.77 at £'o=10 GeV to 
p ~ 0.85 at ^0=10^ GeV. 

In the simulation a correlated pair (oj, (3^) is generated to 



with 



B = 



+p)/2 7(l-p)/2 
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where zi and Z2 are standard normal-distributed random numbers. 

The results of the full and the fast simulations are shown in Fig. 2 and Fig. 3, 
respectively. In these figures the normalized means longitudinal Chcrcnkov 
profiles N{t, At) /A^tot ^-nd their fluctuations (T(r) /N{t, At) for e-induced show- 
ers with different energies are given. A comparison of the SIMEX and FLG 
simulations shows a good agreement in the average profiles. The agreement in 
the fiuctuations in the depth region 0.5 < r < 2 improves with an increase of 
shower energy. 

The average angular distribution of the Cherenkov photons emitted in the 
layer Ar at shower depth r with respect to the shower axis can be written as 
follows: 

P(r, At, cos 9) = N{t, At)'0(t, Ar, cos 9) photons/ster, (19) 
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where 27r j %1){t, At, cos 6) d cos 6—1 
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The function '?/'(r, Ar, cos^) was calculated for each layer Ar = Tj+i —Ti = 0.1 
by the program SIMEX and use for the fast Cherenkov light generation. The 
angular distribution of emitted photons averaged over the cascade depth, 



depends very slowly on the shower energy and is shown in Fig. 4 vs cos 6 for 
the different values of Eq. In this figure the result obtained in Ref.[4] is also 
given for comparison. The angular dependence of iIj{t, At, cos 6) for different 
shower depths is shown in Fig.5. Note that the width of the iIj{t, At, cos 9) 
distribution increases and the average value of cos 9 decreases with r. 



5 Fast generation of the Cherenkov light from shower and hight 
energy muons 

A procedure of the FLGf generation of the Cherenkov photons emitted 
by a shower is the following. The shower is divided into layers with the size 
Ar = 0.1. Each layer is considered as the point-like source of light. Then the 
flux of the Cherenkov photons from the i-th layer at some point X is 



where Ri is the distance between the center of the i-th layer and the point X, 
and 9 is the angle between the shower axis and the direction from the center 
of the i-th layer to the point X. The total light flux at the point X emitted 
by the entire shower is obtained as the sum of contributions from each shower 
layer. 

To check correctness of the simulation using the program FLG, the num- 
ber of photoelectrons collected in the PMT (Hamamatsu R2018) has been 
compared with that from the full MC simulation. The detector conflguration 
simulated for the comparison is given in Fig. (6). The number of the photoelec- 
trons A^pe multiplied by the distance squared between the location of the 
shower maximum and the PMT is shown in Fig. (6) for different orientations 
of the PMT with respect to a shower with energy 1 TeV. The comparison of 
the FLG and the full MC simulation reveals that: 
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• the results are in a good agreement 

• A^pe depends very much on the PMT direction with respect to the shower 

• at the distance R >20 m the shower with energy 1 TeV can be considered 
as a point-hke source of hght. 



The CPU time for the full simulation of the Cherenkov light is large and in- 
creases with the shower energy. But when we use the program FLG the CPU 
time does not depend on the shower energy and decreases dramatically. For 
example, at the energy Eq—IO GeV the CPU time for the FLG simulation 
is 3-10^ times smaller than for the full MC simulation. At the shower energy 
Eo^lO^ GeV the gain is ^S-IO^ times. 

The muon undergoes continuous (ionization) and stochastic energy loss 
processes (5-electrons, bremsstrahlung, e-pair production and muon nuclear 
interaction). In our program for muon tracking the secondary processes be- 
low 0.01 GeV are not simulated individually, but treated as quasi-continuous 
energy loss. The additional light flux A^^ from these low energy processes is 
parametrized by 

N,^No(a + b\nE^). (23) 

Here Nq is the Cherenkov light from the 'naked' muon itself with energy E^. 
The full MC simulation is used for calculation of the Cherenkov light from 
the secondary particles with energies 0.0 1< E <10 GeV. The flux of the 
Cherenkov photons emitted by secondaries with energies >10 GeV is sim- 
ulated by the program FLG. 



6 Conclusions 

A full simulation of secondary electromagnetic processes requires a large 
amount of CPU time, in particular when a low energy threshold for the track- 
ing (f=:ilMeV) has to be used and when large detector volumes are considered. 
In order to speed up the simulation the program FLG has been developed. It 
provides a realistic and fast simulation of the Cherenkov light from electro- 
magnetic showers. Shower to shower fluctuations and correlations are taken 
into account consistently. The agreement between the FLG results and the 
full MC is quite good. The CPU time has been dramatically reduced with 
the FLG simulation. At the shower energy 10^ GeV the gain in CPU time is 
RiS-lO^ times. 

The program for the muon tracking and the FLG codes provide a fast and 
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precise algorithm for large scale Monte Carlo production of high energy events 
for MC study of the underwater/underice Cherenkov detectors like KM3. 
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Figure Captions 



Fig.l: The parameters (a) and of the Gamma distribution vs shower 
energy. The point bars are iTq, and ap. 



Fig. 2: The normahzed mean longitudinal profiles of the Chcrcnkov light N(r, At) /A^tot 
for e-induced showers with the energies Eq = 10, 10^ and 10^ GeV. 



Fig. 3: The relative variance (j{t)/N{t, At) of the Cherenkov light profiles for 
e-induced shower with energies Eq = 10, 10^ and 10^ GeV. 



Fig. 4: The normalized integral angular distribution of the Cherenkov photons 
for different shower energies. 



Fig. 5: The normalized differential angular distribution of the Cherenkov pho- 
tons at different shower depths. 



Fig. 6: * Np,e. as a function of the distance R between the location of the 
shower maximum and the PMT for different directions of the PMT with re- 
spect to the maximum (bold star) at the shower energy E'^IOOO GeV. Arrows 
show the orientations of the PMT: V=cos(6'). 
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Fig.l 
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Fig.3 
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